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1.0  INTRODUCTION 


Current  research  on  flight  control  systems  for  military  aircraft 
emphasizes  digital  control.  Redundant  digital  hardware  is  frequently 
needed  to  achieve  the  required  safety-of-fl ight  reliability,  and  certain 
additional  design  considerations  and  performance  limitations  are  asso- 
ciated with  the  use  of  redundant  hardware. 

When  the  redundant  hardware  is  operated  in  an  asynchronous  manner, 
errors  arise  that  are  attributed  solely  to  asynchronism.  These  errors  are 
present  even  when  the  hardware  is  functioning  in  a normal  manner.  The 
errors  are  small,  but  not  negligible  and  they  influence  dynamic  perfor- 
mance; also,  quantitative  bounds  for  such  "inherent  errors"  are  needed  for 
the  design  of  voting  algorithms  and  to  permit  the  distinction  between  normal 
operation  and  operation  in  the  presence  of  equipment  failure. 

This  report  summarizes  recent  progress  achieved  under  an  Air  Force 
grant  for  the  study  of  inherent  errors  in  asynchronous  digital  flight  con- 
trols. The  research  involves  three  areas: 

1.  Model  Development 

2.  Software  Development 

3.  Parametric  Analysis  of  Representative  Flight  Control  Systems 
The  results  for  each  of  these  areas  are  summarized  below.  A fuller  de- 
scription of  these  results  will  appear  in  an  AFOSR  Scientific  Report,  which 


is  in  preparation. 
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2.0  MODEL  DEVELOPMENT 

Three  models  have  been  developed  for  the  closed-loop  dynamic  operation 
of  an  aircraft  with  dual-redundant  asynchronous  digital  controllers.  All 
three  models  have  the  following  characteristics^:  (1)  Asynchronous  sampling 

is  modeled  by  skewed,  synchronous  sampling  in  which  each  redundant  digital 
controller  operates  at  exactly  the  same  rate  but  there  is  a fixed  time 
skew  for  each  one.  (2)  Although  there  are  two  digital  controllers  only 
one  is  used  to  obtain  the  signals  fed  back  to  control  the  aircraft;  the 
other  controller  does  not  influence  the  aircraft  performance  but  it  is 
used  for  computing  inherent  errors.  (3)  The  aircraft,  sensors,  and  actuators 
are  modeled  by  linear,  continuous-time  state-variable  equations. 

2.1  Model  I 

Figure  1 shows  the  block  diagram  of  the  first  model.  This  model  has  the 
following  set  of  first-order  discrete-time  equations 

x(tk+1)  = F(T,t)  x(tk)  + G(tk,tk+1,x,wp(t)) 

where  x(tk)  is  a combined  state  vector  consisting  of  the  plant  variables  xp 
and  the  digital  controller  variables  xc-j  and  xc2,  as 

yv 

X<V  ■ xcl(tk> 

xc2*y 

F(t,T)  is  a matrix  that  depends  on  T,t,  and  the  aircraft-model  parameters. 


AOC : ANALOG-TO-DIGITAL  CONVERTER 


DAC:  DIGITAL -TO- ANALOG  CONVERTER 

T : SAMPLE  PERIOD 

t:  SKEW 

FIGURE  1 BLOCK  DIAGRAM  FOR  MODEL  I 
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The  inherent  error  for  the  first  model  is  shown  as  e(t)  on  Figure  1. 

This  error  is  piecewise  constant: 

e(t)  = 

The  inherent  error  can  be  characterized  in  a statistical  manner.  Assume 
that  the  external  input  Wp  is  a sample  function  from  a Gaussian  white  noise 
random  process  with  zero  mean.  (In  a typical  flight-control  application  this 
input  could  be  used  to  simulate  the  spectrum  of  frequencies  of  the  actual 
pilot  input  by  first  passing  the  white  noise  signal  through  an  appropriate 
linear  filter;  the  filter  is  chosen  so  that  the  power  spectral  density  of 
its  output  is  representative  of  that  of  a typical  pilot  input;  the  dynamic 
equations  of  the  filter  are  incorporated  into  the  aircraft-model  state 
equations.)  The  covariance  matrix  of  the  states  is  defined  as 

Px(k)  * E[x(tk)  xT(tk» 


' V(tk>  - V>(VT>  for  *k  + T i 1 * Vi 
eB(t)  ' ycl(tk+1>  - for  Vl  i ‘ Vi  + T 


where  E is  the  expected-value  operator.  Px(k)  is  used  in  the  statistical 
characterization  of  the  inherent  errors. 

The  steady-state  covariance  of  the  states,  pxss»  is  found  by  solving 
the  equation 


Pxss  = F(T’t)  Pxss  F(T*t)  + V(T’t) 

where  V(T,x)  Is  a matrix  that  can  be  calculated  from  system  matrices  and  the 
covariance  matrix  for  Wp.  Once  Px$$  is  obtained,  the  steady- state  covariance 


matrices  of  the  inherent  errors,  e^  and  eg,  can  be  obtained  by  matrix 
multiplication.  The  numerical  values  of  the  elements  of  Px$$  give  the 
variances  of  the  inherent  errors  about  their  means  of  zero. 

All  required  equations  for  the  state-variable  model  and  the  statisti- 
cal analysis  have  been  developed  and  are  given  in  References  1 and  2. 

2.3  Model  II 

The  second  model  incorporates  the  same  assumptions  as  the  first  model 
but  also  allows  for  three  time  delays,  associated  with  the  delays  due  to 
data  transmission  and  control-signal  computations.  Model  II  assumes  sampling 
at  a single  rate  and  includes  sampling  of  the  pilot  input. 

The  three  time  delays  are: 

6,:  A/D  conversion  and  transfer  time 

d 

<$c:  digital-controller  computation  time 

/ 

6^:  D/A  conversion  and  transfer  time 

Figure  2 shows  the  basic  configuration  assumed  and  Table  1 gives  the  event 
times  where  relevant  changes  in  the  system  variables  occur. 

The  state  equations  for  Model  II  involve  the  new  variables  and  x^, 
which  are  needed  to  provide  the  delay  in  the  piecewise-constant  samples  of 
the  Xp,  for  digital  controllers  1 and  2,  respectively.  Thus  the  state 
vector  for  Model  II  is 


ADC:  ANALOG -TO- DIGITAL  CONVERTER 

DAC  : DIGITAL-TO- ANALOG  CONVERTER 

S/H  : SAMPLE-AND- HOLD 


T : SAMPLE  PERIOD 

T : SKEW 


FIGURE  2 BLOCK  DIAGRAM  FOR  MODEL  XL 


Table  1.  Event  Times  for  Model  II 


Events 


+ 6a  + 5c 


yp  sampling  initiated  for  digital  controller  1 
Wp  sampling  initiated. 

yp  reaches  digital  controller  1. 
yc1  computation  is  complete. 


+ 6a  + 5C  ♦ Jd  *cl 


reaches  the  actuator. 


+ x + 6. 


yp  sampling  initiated  for  controller  2. 
yp  reaches  controller  2. 


+ t + + 6C  yc^  computation  is  complete. 
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and  the  state  equations  are 


x(tk+1)  = F(T,t)  x(tk)  + G(T,t)  w (tk) 


The  matrices  F and  G in  this  equation  are  different  from  those  in  the  state 
equation  for  Model  I;  in  particular,  G(T,x)  is  dependent  on  T and  x but  not 
on  wp(t)  since  the  pilot  input  is  sampled  in  Model  II. 

The  inherent  error  is  again  piecewise-constant  and  is  defined  as 


e(t)  = 


e^(t)  = ycl(W6c)  " yc2^tk+T+6a+6c^  for  t in  Interval  A 


(eB(t)  . yci  (tk+i+{a+5c)  - yc2(tk+T+lSa+<Sc^  ^or  ^ 1n  Nerval  B 


where  Interval  A is  [tk+x+6a+6c,  ^k+l+5a+l5c^  anc*  Interval  B is 

^•tk+l+<Sa+<5c’  tk+l+T+6a+6c^ ' k = • 

The  covariance  analysis  for  Model  II  follows  the  same  procedure  that  was 
used  in  Model  I.  The  covariance  matrix  for  the  states  is  defined  as 

Px(k)  * E[x(tk)  xT(tk)] 

and  it  satisfies  the  difference  equation 

P (k+1)  = FP  ( k ) FT  + GWGT 

A A 

The  covariance  matrices  for  e^(t)  and  eB(t)  are  related  to  Px(k)  as  follows 
PeA(k)  a (Hr»2)  Px(k)  (HrH2)T  + pWpT 
PeB(k)  ° (HiF-H2)  Px(k)  (H-|F-H2)T  + (HjG-p)  W^G-p)1 
where  Hjt  H2,  and  p are  constant  matrices. 


2.3  Model  III 

The  third  model  incorporates  all  the  assumptions  of  the  first  model 
but  also  allows  for  multirate  sampling;  that  is,  the  pilot  input  wp  is 
now  sampled  at  a slower  rate  than  the  rate  at  which  the  control  system 
operates.  Figure  3 shows  the  configuration  for  Model  III.  Note  that  the 
digital  controller  sampling  rates  are  n times  the  pilot-input  sampling 
rate,  (n  is  a positive  integer.) 

The  state  equations  for  Model  III  take  the  form 


x(tk+1)  = F(T,t)  x(tk)  + G(T ,t ) wp(tk) 


where  F and  G are  different  from,  but  play  the  same  role  as  F and  G in  the 
equations  for  Models  I and  II. 

The  inherent  error  is  again  piecewise-constant,  and  it  changes  2n  times 
in  the  time  interval  from  tk+x  to  tk+T : 


e(t)  = 


- J'c2(tk+Wr*T)  for 
' ycl(tk*TLT)  - yc2<tk4T+t>  f°r 


i = 0,1 n-1,  k = 0,1,2,...  . 


It  is  possible  in  this  case  to  define  "average  errors"  and  Eg,  as 


, n-i 

EA  = n E eAi 


Ee  ‘ n ef 
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ADC:  ANALOG -TO-  DIGITAL  CONVERTER 

DAC  : DIGITAL -TO-  ANALOG  CONVERTER 

S/H  ! SAMPLE- AND- HOLD 

T : PILOT -INPUT  SAMPLE  PERIOD 

T/n  : DIGITAL  - CONTROLLER  SAMPLE  PERIOD 

T : SKEW 


FIGURE  3 BLOCK  DIAGRAM  FOR  MODEL  HI 
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and  to  obtain  closed-form  expressions  for  both,  in  terms  of  x(tk)  and 

vv- 

The  covariance  analysis  for  Model  III  is  complicated  by  the  presence  of 
multiple  rates.  Define  the  covariance  matrix  of  the  state  variables  by 

Px(k,m)  = E[x(tk  + £T)  xT(tk  + JT)] 
for 

m = 0,1,...,  m-1 , k = 0,1,2,...  . 

Then  it  can  be  shown  that 

Px(k,m+1)  = AF  Px(k,m)AFT  + AGW^G1 

where  AF  and  AG  are  known  functions  of  the  system  parameters  and  sample 
periods  and  W is  the  input  covariance  matrix.  The  above  equation  requires 
a starting  value,  Px (k ,0 ) , which  may  be  obtained  from 

Px(k+1,0)  = F(T,t)  Px(k,0)FT(T,x)  + G(T,t)U0GT(T,t) 

These  equations  are  also  useful  for  studying  the  steady-state  covariance 
matrices  for  the  state  variables. 

The  2n  covariance  matrices  for  the  inherent  errors  are  related  to  the 
covariance  matrices  for  the  states  in  a straightforward  manner. 


3.0  SOFTWARE  DEVELOPMENT 


The  use  of  the  digital  computer  is  needed  to  use  the  models  for  dynamic 

simulations  and  covariance  analyses  of  realistic  flight-control  systems.  A 

few  low-order,  simplified  closed-loop  systems  have  been  analyzed  by  hand 
1 2 

with  Model  I,  * but  practical  systems  must  be  studied  by  using  specially 
written  software  packages.  There  is  a separate  package  for  each  model. 

These  packages  are  described  briefly  below.  Details  for  the  package 
for  Model  I are  in  References  1,  and  are  in  preparation  for  Models  II  and  III. 

3. 1 Model  I Software 

The  software  package  for  Model  I consists  of  a FORTRAN  main  program  and 
1 2 

16  subroutines.  ’ There  are  approximately  1500  cards  in  the  package.  Most 
of  the  subroutines  are  taken  from  the  software  package  DIGIKON,  which  was 
written  under  contract  with  the  Air  Force  Flight  Dynamics  Laboratory  by  the 

3 

Honeywell  Flight  Systems  Division. 

The  basic  input  to  the  software  package  is  the  state  variable  description 
of  the  aircraft  dynamics  and  the  digital  controllers.  Output  consists  of  the 
description  of  the  closed-loop,  discrete-time  state  equations,  the  steady- 
state  covariance  matrices  for  the  states  and  the  inherent  error,  and  the 
system  response  to  a unit-step  pilot  input. 

3.2  Model  II  and  Model  III  Software 

There  are  distinct  software  packages  for  Model  II  and  for  Model  III. 

They  are  similar  in  structure  and  function  to  the  software  for  Model  I and 
make  use  of  the  same  basic  DIGIKON  subroutines. 


>* 
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4.0  PARAMETRIC  ANALYSES 

The  practical  system  chosen  for  study  using  the  models  and  software 

previously  described  is  the  pitch  axis  flight-control  system  of  the  A-7D 

tactical  fighter.  This  aircraft  is  the  one  whose  dynamics  are  being 

simulated  in  the  Flight  Control  Facility  of  the  Air  Force  Flight  Dynamics 

Laboratory's  Digital  Avionics  Information  System  (DAIS)  Advanced  Develop- 
4 

ment  Program.  The  purpose  of  the  DAIS  program  is  to  demonstrate  a co- 
herent solution  to  the  costly  problem  of  proliferation  and  nonstandardiza- 
tion of  aircraft  avionics  through  modular  digital  hardware  and  software 
concepts.  The  program,  managed  by  the  Air  Force  Avionics  Laboratory,  is 
a joint  effort  of  several  Air  Force  organizations  and  includes  the  Air 

- t 

Force  Flight  Dynamics  Laboratory. 

References  1 and  2 contain  discussions  of  the  results  of  applying 
Model  I to  a version  of  the  A-7D  pitch  axis  flight-control  system.  The 
major  system  components  include  the  pilot-input,  the  aircraft  and  actuator 
dynamics,  and  dual -redundant  asynchronous  digital  controllers.  The  first 
two  of  these  components  comprise  a 4th-order  linear,  continuous  system  and 
the  controllers  are  each  first-order. 

The  results  of  the  study  indicate  that,  under  typical  conditions  of 
flight  (nominal  forward  velocity  in  the  cruise  condition  at  sea  level)  the 
absolute  value  of  the  largest  inherent  error  in  the  elevator-deflection 
command  is  3.20  percent  of  the  steady-state  value.  When  a zero-mean, 
white-noise  input  is  applied,  the  covariance  analysis  revealed  a steady- 
state  inherent  comparison  error  of  3.13  percent  maximum.  See  References  1 
.and  2 for  a more  complete  discussion. 
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Work  Is  progressing  on  applying  all  three  models  to  a slightly  modified 

version  of  the  A-7D  pitch  axis  control  system.  This  version  follows  the 

equations  programmed  in  analog-computer  simulation  that  is  currently 

5 

being  used  in  the  DAIS  program.  The  linear  dynamics  are  6th  order  and  the 
controller  equations  are  again  first-order.  These  studies  will  be  reported 
in  the  AFOSR  Scientific  Report  which  was  mentioned  in  Section  1.0  and  which 
is  in  preparation. 


' --  - -- \ • 


16 


5.0  REFERENCES 


1.  V.J.  Darcy  and  C.  Slivinsky,  "Analysis  of  Inherent  Errors  in  Asynchronous 
Redundant  Flight  Control  Systems,"  Air  Force  Flight  Dynamics  Lab.,  Air 
Force  Systems  Command,  Wright-Patterson  Air  Force  Base,  OH,  Tech.  Rep. 
AFFDL-TR- 76-16,  1976. 

2.  V.J.  Darcy  and  C.  Slivinsky,  "Inherent  Errors  in  Asynchronous,  Redundant 
Digital  Flight  Controls,"  IEEE  Transactions  on  Automatic  Control, 

February,  1977. 

3.  A.F.  Konar,  et  al_. , "Digital  Flight  Control  Systems  for  Tactical  Fighters, 

Vol.  I:  Digital  Flight  Control  Systems  Analysis,  Vol.  II:  Documentation 

of  the  Digital  Control  Analysis  Software  (DIGIKON),"  Air  Force  Flight 
Dynamics  Laboratory,  Air  Force  Systems  Command,  Wright-Patterson  Air  Force 
Base,  OH,  Tech.  Rep.  AFFDL-TR- 75  119,  Vol.  I and  II,  June,  1974. 


L 


4.  AFFDL/DAIS  Flight  Engineering  Faci 1 Program  Plan  prepared  by  TRW 
Defense  and  Space  Systems  Group,  for  Air  Force  Flight  Dynamics 
Laboratory,  Air  Force  Systems  Command.  ight-Patterson  Air  Force  Base, 
Ohio,  April,  1977. 

5.  Martin , J.M. , "Analog  Computer  Simulation  of  the  A-7D  Aerodynamic  and 
FCS  for  DAIS,"  Air  Force  Flight  Dynamics  Laboratory,  Air  Force  Systems 
Command,  Wright-Patterson  Air  Force  Base,  Ohio.  To  appear. 


i ’ ~ y ",~;>T~r-  "*' 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (When  hntu  h.nterod) 


IEPORT  DOCUMENTATION  PAGE 


1.  REPORT  NU^feC 

AFOSR&R- 


jif>i:  READ  INSTRUCTIONS 

OL~ HEI  ORK  COMPLETING  FORM 

2.  GOVT  ACCE'  -’IOM  NO.  3-  RECIPIENT'S  CATALOG  NUMBER 


ANALYSIS  OF  INHERENT  ERRORS  IN 
^DIGITAL  FLIGHT  CONTROLS  . 


ASYNCHRONOUS  ( 


1 7.  author^*; 


Charles  /Slivinsky 


@JTd  'Y 


Interim 


PIING/>^G.  REfORT  NUMBER 

-£1 


AFOSR  76-2968 


: 

j 


9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  S*  \ 

University  of  Missouri  QiJ- — -< 

Department  of  Electrical  Engineering  j V rt  1 
Columbia,  Missouri  65201 f~  P I 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB,  Washington,  DC  20332 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

A'Wfl'b5?'RK  UNIT  NUMBERS 


'A1  0Vf\ 


12.  REPORT  DATE 

1977 


Bolling  At  B,  Washington,  DU  ia.  number  of  pages 

16 

1«  MONITORING  AGENCY  NAME  » ADDRESSf/f  different  from  Controlling  Otlica)  15.  SECURITY  CLASS,  (ol  this  rrport) 

I 7”: - UNCLASSIFIED 


r/r  H F 


1 5a.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  <hfs  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


1 19.  KEY  WORDS  ( Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 


ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 


Progress  has  been  achieved  in  the  study  of  inherent  errors  in  asynchronous 
digital  flight  controls.  The  research  involves  model  development, 
software  development,  and  parametric  analysis  of  representative  flight  control 
systems . 

\ 


DD  1 JAN  73  1473  EDITION  OF  1 NOV  SS  IS  OBSOLETE 


UNCLASSIFIED  % / / 

SECURITY  CLASSIFICATION  OF  THIS  B AO l fWhw.  DM.  EntarMl)  . 


